Decay by tunneling of Bosonic and Fermionic Tonks-Girardeau Gases 
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We study the tunneling dynamics of bosonic and fermionic Tonks-Girardeau gases from a hard 
wall trap, in which one of the walls is substituted by a delta potential. Using the Fermi-Bose map, 
the decay of the probability to remain in the trap is studied as a function of both the number of 
particles and the intensity of the end-cap delta laser. The fermionic gas is shown to be a good 
candidate to study deviations of the non-exponential decay of the single-particle type, whereas for 
the bosonic case a novel regime of non-exponential decay appears due to the contributions of different 
resonances of the trap. 
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I. INTRODUCTION 

Decay of a metastable system via tunneling is one of 
the most remarkable and old effects in quantum mechan- 
ics. Since Gamow's analysis of alpha decay, resonance 
theory, which applies to virtually all fields from particle 
to molecular physics, has been motivated by this phe- 
nomenon. Simple treatments examine the escape or sur- 
vival of single particle wave functions in one dimensional 
(ID) potentials. At this level much attention has been 
paid to deviations from exponential decay, and Zeno or 
anti-Zeno effects. Also, exact results are typically avail- 
able, by means of analytical models or numerically. The 
more complex decay of a multiparticle unstable system is 
treated by more sophisticated multichannel, or reactive- 
scattering approaches, sometimes with statistical approx- 
imations or, depending on the system and environment, 
in a phenomenological way, and also using mean-field ap- 
proximations. In "macroscopic quantum tunneling" , a 
macroscopic variable, such as the phase difference of the 
Cooper pair wave function across a Josephson junction 
obeys a simple tunneling equation for an effective par- 
ticle subjected to dissipation 0,0- The effect of dissi- 
pation due to the perturbing environment has thus been 
extensively discussed and measured. Other macroscopic 
quantum tunneling effect much studied in recent times 
is the tunneling and decay of Bose-Einstein condensates; 
in particular the effect of effective atom-atom interaction 
and the non-linear term in the Gross-Pitaevskii mean- 
field approach HSU- 



'Electronic address: qfbdeeca@ehu.es 
^Electronic address: qfbdeacf@ehu.es 

t Permanent address: Instituto de Fisica, Universidad Nacional 
Autonoma de Mexico, Apartado Postal 20 364, 01000 Mexico, D.F., 
Mexico; Electronic address: gaston@fisica.unam.mx 
§ Electronic address: jg.muga@ehu.es 
^Electronic address: raizen@physics.utexas.edu 



Some works go beyond the mean field theory using 
simplified Hamiltonians 6], which is particularly rele- 
vant for few-body systems. These systems with not too 
many particles may still be amenable of exact treatments 
but will show differences from single and many-particle 
ones. The experimental study of few-body tunneling is 
a formidable prospect, as it requires initial preparation 
of a ground few-body Fock state, precise control over 
the tunneling time, and the ability to count single atoms 
with unit quantum efficiency. Until recently, such ca- 
pabilities did not exist so that any experimental tests 
seemed unlikely. However recent developments now open 
the door for few body tunneling experiments and moti- 
vate the present work in anticipation of such results. The 
starting point was the development of a novel optical box 
trap that confines a degenerate Bose gas, together with 
single-atom counting Eg . The same box trap was used 
to produce number squeezing of atoms by confining a de- 
generate Bose gas and controlled lowering of the walls un- 
til a final value @ . The observed fluctuations in number 
were a factor of two below the Poissonian limit, but the 
residual noise can be accounted for by known sources of 
technical noise, so that these experiments are consistent 
with number-state production. This simple procedure, 
called "laser culling of atoms" has recently been analyzed 
theoretically, and is shown to produce atomic few-body 
Fock states for sufficiently slow ramp time and neglecting 
quantum tunneling through the barrier The latter 
effect can be highly suppressed by sculpting the shape of 
the barrier using the techniques described in 8] . The bar- 
rier width can then be reduced at a well defined time, al- 
lowing quantum tunneling to occur. This system should 
therefore enable the first experimental study of few body 
quantum tunneling in different regimes of interaction and 
with controlled number. 

In this paper we study the decay from a trap by tunnel- 
ing trough a delta barrier of a few-body Tonks-Girardeau 
gas with the aim of obtaining exact results of a few-body 
decay problem for which the experimental verification is 
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m view m,|£J,|2|. In particular, we find few-body devia- 
tions from the exponential decay law. For one particle, 
deviations from exponential decay have been long pre- 
dicted at both short and long times. Short time devia- 
tions were observed experimentally with ultracold atoms 
[Tlj . whereas long time deviations have been observed 
very recently for the first time in dissolved organic ma- 
terials |l~2]. A suitable system to test deviations from 
exponential decay at the few-body level is the bosonic 
Tonks-Girardeau (BTG) gas E3- Such a gas actually 
mimics the fermionic behavior to minimize the strongly 
repulsive interaction. At low densities, the TG regime 
can be reached under a strong enough radial confinement 
[T3 | such that the transverse degrees of freedom are re- 
duced to zero-point oscillations, resulting a ID effective 
system. Indeed, experimental observations have already 
been reported [lq[. 

As another striking example of the Fermi-Bose duality 
[H |. the fermionic Tonks-Girardeau (FTG) gas has been 
described as a ID spin-aligned fermionic gas in the limit 
of highly attractive interactions mediated by a 3D p-wave 
Feshbach resonance [l^ ■ The FTG gas is in a sense the 
opposite of the BTG gas and exhibits "bosonization" be- 
cause of the strong attraction. 

In this work we study the tunneling dynamics of both 
bosonic and fermionic TG gases, initially confined in a 
hard wall potential. For the bosonic case, the Fermi-Bose 
(FB) map El El El El relates the wavefunction of N 
strongly interacting bosons, tps, to the one of an ideal 
Fermi gas with all spins frozen in the same direction, ip F , 
where the tilde indicates the "dual" , or auxiliary system. 
The Fermi wavefunction, is built as a Slater determinant, 
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where (j) n (xk) are the eigenstates of the trap. To account 
for the proper quantum statistics, the "antisymmetric 
unit function" 



A = Y\ sgn(xk — Xj) 

l<j<k<N 

is introduced in such a way that 



(2) 



ip B (xi,...,x N ) = A(x%, . . . ,x N )ip F (xi, . . . ,x N ). (3) 

In so doing, iPb(xi, ■■■,%n) becomes totally symmetric 
under permutation of particles. Moreover, since A is in- 
volutive and does not include time explicitly, for any uni- 
tary evolution it can be proved that the time-dependent 
density profile can be calculated as Ell 



p B (x,t) = N \ip B {x,X2, ■ ■ ■ 1 x N \t)\ 2 dx 2 



■ dx 
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In a similar way, one can deal with the Fermionic TG 
(FTG) gas E3." 

but using now the generalized FB map- 
ping in the opposite direction. The wavefunction of the 
FTG gas can be written in terms of the Hartree product 
describing the dual (auxiliary) system, which is now the 
ideal Bose gas, 
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iPb(xi,- ■ - ,x N ) = Y[4>o(xi). 



(•5) 



The proper symmetrization is carried out also through 
the antisymmetric unit function, 

ip F (xi, - ■ ■ ,x N ) = A(xt,..., x N )^ B {xi,...,x N ), (6) 

in such a way that the density profile can be written as 

p F (x,t) =N\Mx,t)\ 2 . (7) 

The trap we will consider is of hard-wall kind after its 
experimental realization in its all-optical 20] or micro- 
electronic chip version |2l|. These traps have recently 
generated much theoretical work in the field of ultracold 
atoms in low dimensions El El El El . 



II. SINGLE EIGENMODE DYNAMICS 

First of all we shall study the time evolution of the 
n-th eigenstate of a hard wall trap. As it is well-known 
they have the general form 



/2 / T17TX \ 

^ sin y—jr-J X[o,L](x), 



(8) 



with n 6 N and X[o l]{x) the characteristic function in 
[0,L\. 

At time equal zero the right-wall is substituted by a 
delta potential, V(x) — r)'5(x — L), which is to repre- 
sent a far-detuned laser from the atomic resonance, so 
that atomic excitation becomes negligible. This model 
has been recently considered to study non-exponential 
decay at both short and long times at the single-particle 
level in |2^. The time evolution of a given initial state 
may be written in terms of the retarded Green's function 
g(x, x'; t) as 



(j) n (x,t) = / g(x,x';t)(j) n (x',0)dx'. 



(9) 



Equation @ may be calculated using an expansion in 
the eigenfunctions of the Hamiltonian {|fc + )|fc G 
i.e., the so called physical wave solutions, 

/>oo 

<f> n (x,t) = / dkix^ik+l^e-^ 2 ^, (10) 
Jo 

r— ( sin(fcr)/ J+(k), x < L 
(x\k+) = J-\ (11) 
V n [ (i/2) [e~ lkx - S{k)e lkx ] , x > L, 
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FIG. 1: Diagram on the complex fc-plane of the first five 
resonance and anti-resonance poles of the hard-wall leaking 
trap at two different strengths of the end-cap laser. Here 
and in all figures we use dimensionless units with L = 1 and 
2m = h = 1. 



FIG. 2: Single particle nonescape probability exhibiting the 
long time nonexponential behavior for the first ten eigenstates 
of the hard wall trap, tunneling through a delta potential with 
r\ = 5. From top to bottom: n = 1, . . . , 10. 



where the S-matrix S(k) = J_(fc)/J + (fc), with J-(k) = 
J+(k), and the Jost function J+(k) reads 

J + (k) = l + ^L(a 2ikL -l). (12) 

Above and in the rest of this work we take k = 
[2mE /ft 2 ] 1 / 2 and r] = [2m/ft 2 ]r?', E being the energy of 
the decaying particle and m its mass. 

The position of the pole kj of the S'-matrix in the k- 
plane gives us information about the lifetime h/Tj of the 
resonance and its real energy £j in the leaking trap (see 
Fig. 1): Ej = k 2 ] h 2 /{2m) = Sj - iTj/2. We notice that 
for a weaker laser (smaller 77), the resonance decays faster. 
Besides, j is chosen so that the width of the resonance 
increases monotonically with it, j = 1 corresponding to 
the longest lived resonance. 

The limit rj — > is an important reference correspond- 
ing to free evolution in the presence of a wall at the origin, 
once the endcap laser at x — L has been turned off. In 
such fully analytical solution is available invoking 

the method of images [22j], see a detailed study of the 
BTG expansion dynamics in the absence of the wall in 

The decay of the different eigenstates of the hard wall 
trap is shown in Fig. [3 It exhibits the characteristic 
transition to long times as an inverse power of time. The 
scale of the figure conceals an interesting behavior at 
short times, which is exhibited in Fig. [3] We observe, 
with exception of the first decaying state, that each state 
presents some characteristic oscillations at short times 
and then tends to decay with the same slope as the first 
state. Physical insight regarding this behavior may be 
obtained making use of the formalism of resonant states. 
In this representation one Laplace transforms g(x,x';t) 
in Eq. @ into the complex fc-plane to exploit the ana- 



lytical properties of the corresponding outgoing Green's 
function G + (x, x'\ k) j^. This leads to an alternative 
expression to Eq. I|ll|l . namely, 

00 
3=1 

- ( cj} n (x',0)dx' ( G + (x,x';k)e- ihk2t ^ m kdk, 

n JQ Jc l 

(13) 

where the sum runs over the so called proper complex 
poles i.e., those on the fourth quadrant of the complex 
fc-plane with Re(fcj) > |lm(fcj)|; the Uj's obey outgoing 
boundary conditions at x — L and satisfy the Schrodinger 
equation with complex eigenvalues Ej = £j — iTj/2; 
the coefficients Cj (n) give the overlap of the initial state 
(j) n with the resonant states uj of the problem, namely, 

Cj(n) — 4> n (x, 0)uj(x)dx; the integral term involving 
G + (x, x'; k) stands for the non-exponential contribution, 
which in general may be neglected except at ultrashort 
or very long times. The path Cl of the integral is chosen, 
without loss of generality, as a straight line 45 degrees off 
the real fc-axis along the complex fc-plane passing through 
the origin k = [28|. 

The nonescape probability is defined as 

P»(t) = / \Mx,t)\ 2 dx, (14) 
Jo 

(to be distinguished from the "survival probability" 
\{<t>n{t)\<t>n{fy)\ 2 1 which is harder to measure). Using the 
exponential contribution to Eq. (|13[) into Eq. I|14[) yields 
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FIG. 3: Single particle nonescape probability for the first ten 
eigenstates of the hard wall trap, tunneling through a delta 
potential with 77 = 5 at short times. From top to bottom: 
n = 1,...,10. 



the nonescape probability for each decaying eigenstate as 

N 



(15) 

where Ij 8 — J Uj(x)u*(x)dx. The exponential term in 
Eq. I|15l) . with N — n, reproduces very well the regime 
depicted in Fig. [3] It reflects a transient regime where 
each eigenstate 4> n makes eventually a transition into the 
longest lived eigenstate (j — 1). As shown in Fig. 
after the exponential regime, it follows the long time t~ 3 
inverse power law governed by the integral contribution 
to Eq. (|13|l . Standard asymptotic analysis gives the result 



plon 9{t) 



2m\ 3 L 3 C 2 (n 
~YJ 12Tr{l + T 1 L) 4 ~t 3 ~ 



(16) 



where C{n) = J c ^ <f> n (x,0)xdx = LV2L(-l) n /(mr). 
Clearly the transients in the transition from exponential 
to nonexponential behavior observed in Fig. [21 originate 
from the interference between the exponential and long 
time expressions of 4> n {x,t). The transition may be dis- 
placed to earlier times and made more easily observable 
by decreasing the ratio R = £\/T\ of the longest lived 
resonance [M, 301. 



III. TONKS-GIRARDEAU GAS 

We next generalize the notion of nonescape probability 
to a few particle system, associating it with the average 
number of particles within the trap at a given time, 
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FIG. 4: Logarithm of the average number of particles within 
the trap for a BTG gas as a function of time, for different 
values of the total number of particles. From top to bottom, 
N = 1, . . . , 10 with 77 = 5. The curve for the FTG gas of 
iV-particles is obtained by shifting up the N = 1 BTG gas 
curve (single-particle case) by lniV. 



where p{x, t) is the density normalized to the total num- 
ber of particles N . The nonescape probability per par- 
ticle thus becomes P(t) = Nt/N. From this definition 
it is clear that for the FTG gas, P(t) is identical to the 
single particle nonescape probability associated with the 
ground state of the trap. However, this means that the 
total signal is enhanced by a factor corresponding to the 
number of particles, a key advantage for the experimen- 
tal study of deviations from the exponential decay law. 
In other words, In Nt for an iV-particle FTG gas in its 
ground state, is obtained by shifting upwards the curve 
for n = 1 in Figs. Hornby ^N. 

More remarkable yet, the BTG gas exhibits novel few- 
body decay features. In particular, a new regime of non- 
exponential decay arises from the sum of contributions 
of the single particle exponentials. Figure 0] shows that, 
for the BTG gas, the higher the number of particles the 
faster is the tunneling rate and, as a result of the differ- 
ent energy contributions, the short time dependence is 
affected. 

The effect of the end-cap laser intensity can be simu- 
lated by varying 77. In so doing, we have learnt that the 
resonance poles shift towards the real axis in the fc-plane 
for increasing 77, see Fig. ^ For vanishing 77, even for the 
single particle case, and therefore the FTG gas, there is 
no reason to expect exponential decay since the potential 
is not a trap any more and there are no resonances. 

This behavior is shown in Fig. [5] For a given 77, a 
clear difference between BTG and FTG gases is the decay 
rate, which is larger for the former at short times, a direct 
consequence of the different ground state energy for both 
systems, namely, 



and 



N T (t) = / dxp(x,t), 



(17) 



Ebtg = h 2 n 2 N(N + l)(2N + l)/(12mL 2 ) 



E FTG = h 2 n 2 N/{2mL 
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FIG. 5: Logarithm of the average number of particles in the 
trap for a BTG gas (dashed lines) and FTG gas (continuous 
lines) of 10 particles for different intensities of the end-cap 
laser (for both cases, from top to bottom 77 = 10, 5,0). 

However, the high resonance contributions eventually di- 
minish, and the rates become finally equal being governed 
by the longest lived resonance, deep into the purely (sin- 
gle) exponential regime. Notice that at short times, the 
decay is dominated by high energy components and it 
is essentially independent of the strength of the delta 
potential. Moreover, at short times the non-exponential 
deviation is concave-up for the BTG system and concave- 
down for the fcrmionic case. 

Regarding the long time behavior, it follows from Eq. 
(JTS) that N T (t) = NP (t) for the FTG gas exhibits the 
1/t 3 dependence. This is also the case for BTG gas where 



tf r (t) = $>»w«£T- (18) 

n— 1 n— 1 

This result guarantees that the generalization for our few- 
body systems of the dwell time (the average time spent 
by a particle in a spatial region) is possible, being a mean- 
ingful finite quantity ■ 

In conclusion, we have studied exactly two related few- 
body tunneling problems, namely, that of bosonic and 
fermionic TG gases. The FTG gas has been pointed out 
as a good system to observe long time deviations from 
exponential decay because of the strength of the signal, 
proportional to the number of particles. The recent first 
measurement of long time deviations [l2 | is based on the 
important effect of the environment in organic molecules 



in solution, so that the deviation for "pure" , isolated sys- 
tems remains to be observed [3(3 ■ For the bosonic case 
a new deviation of the exponential decay appears, which 
can be understood as a sum of N single-particle contri- 
butions. 

It is still an experimental challenge to get to the strong 
TG limit in a flat box, the main limiting factor being the 
confinement in transverse directions. Assuming, accord- 
ing to Reichel and Thywissen |33j| . a "maximal practical 
value" of transversal frequency of 1 MHz, N = 10, a box 
of 10/im, and the constants for rubidium 87 (scattering 
lenght a m 5 nm), the ratio a of interaction energy to the 
potential energy is a ~ 440, whereas the ratio of chemical 
potential to the kinetic energy is 7 « 88, see [33J. For ru- 
bidium 85 at a Feshbach resonance the scattering length, 
a and 7 may increase by a factor of 100 |34| , making op- 
tical traps a viable alternative to magnetic confinement. 
These values point out at a TG gas regime not too far 
away from current capabilities, but the actual implemen- 
tation may still be difficult also because of the need for 
accurate single atom detection. 

We may in any case expect that reaching a strict TG 
regime is not absolutely essential to find interesting few- 
body decay effects. The upshot of our analysis is that 
for the few-body ground state of the box, the BTG gas 
exhibits maximal deviations from the exponential law at 
short times, whereas they are completely absent in the 
FTG gas or ideal Bose gas. The behaviour between the 
ideal Bose gas and BTG gas could be smoothly extrapo- 
lated by increasing the effective interactions [3J|. More- 
over, the Bose- Fermi duality pl| opens up the possibility 
of observing experimentally such effects in a wider class 
of systems. 
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